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In this work a method is presented to derive the generating functional in path 
integral form for a system with an arbitrary number of degrees of freedom and 
constrained by general conditions. The method is applied to the case of the dynamics 
of an inextensible chain subjected to external forces. Next, the generating functional 
of the inextensible chain is computed assuming that the interactions are switched 
off. Finally, the generating functional of a two dimensional nonlinear sigma model 
with 0(3) symmetry is derived exploiting its similarities with the model describing 
the dynamics of the inextensible chain. 

INTRODUCTION 

Stochastic systems in the presence of constraints are studied in connection with several 



physical problems [l-|9|. In the absence of constraints, a very useful method in order to 



solve stochastic differential equations 
their solution in path integral form 



ike those of Fokker-Planck consists in expressing 



10l4l3l|. Once this has been done, a rich choice of path 
integral techniques is available for concrete calculations. The derivation of the path integral 
solutions of the Fokker-Planck equations becomes however complicated when constraints are 
added. Various approaches have been proposed in order to cope with constrained stochastic 



differential equations j8l.l9l.fi4 



30| . but the available path integral formulations for statistical 



systems out of equilibrium still remain awkward to be used in practical applications. 

In this work a method is presented to derive for a system out of equilibrium and with 
an arbitrary number of degrees of freedom a path integral expression of the generating 
functional of the relevant correlation functions. The starting point are the Langevin equa- 
tions describing the unconstrained system with the addition general constraint conditions. 
Roughly speaking, the strategy behind the method is to introduce fictitious degrees of free- 
dom and white noises in order to allow the variables that should be constrained to fluctuate. 
It turns out that the dynamics of this extended system can be described by an enlarged set 
of overdamped Langevin equations, so that it is possible to construct the path integral of the 
related generating functional by using standard techniques. The original system of Langevin 
equations and of constraints is recovered in the limit in which the friction coefficients of the 
fictitious degrees of freedom and the widths of the fictitious noise distributions approach 
zero. In this limit, in the path integral of the generating functional appear functional Dirac 
delta functions whose role is to impose the constraints. The proposed method is tested in 
the case of a discrete inextensible chain, for which the generating functional of the rele- 



vant correlation functions has already been derived in path integral form in 
undergone several consistency checks 



32 



3l| and 



has 



33| . A considerable advance with respect to 



3l| 



is the inclusion in the treatment of the external forces, that were previously missing. While 
with the usual approaches the path integral of t he g enerating functional contains extra de- 



8, 



9, 
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35| or, alternatively, is complicated by 



grees of freedom (the Lagrange multipliers' 
functional determinants of block matrices 8|, 19|, 136| , in the present approach the functional 
determinants are relatively simple and can be computed analytically. 

Even with the simplest possible formulation, the treatment of system with constraints is 
complicated by the fact that constraints involve nonlinear functions which do not contain 
small or large parameters. In this situation, any approximation implies a change of the 
imposed conditions. For instance, within the present approach the main problem are the 
functional delta functions imposing the constraints. They make the computation of the gen- 
erating functional extremely complicated. Usually, the constrained path integral is reduced 
to a gaussian one at the price of drastic simplifications 
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. In the second part of this 



article it is shown in the case of the inextensible chain that it is possible to go beyond the 
gaussian approximation. To this purpose, the path integral of the generating functional of 



3 



an inextensible chain is rewritten as a product of two path integrals describing two systems 
interacting with an auxihary field. One of the two path integral is gaussian, while the second, 
which contains the Dirac delta function coming from the constraints, may be computed by 
performing a series expansion in powers of the auxiliary field. The series converges provided 
a parameter which is the analog of the Planck constant in quantum mechanics is small. An 
exact expression of the terms of any order appearing in the series is derived. The generating 
functional of these terms is computed in closed form. 

The path integral of a inextensible chain in three dimensions is very similar to that of 
an 0(3) nonlinear sigma model 40|] in two dimensions. For this reason, it is possible to 
apply also to this type of sigma models the same method used to evaluate the generating 
functional of the inextensible chain. 

The material presented in this paper is divided as follows. In the next Section, a method 
for constructing the generating functional of a general constrained stochastic system is pre- 
sented. As a concrete example, we apply this method to the case of a inextensible chain 
consisting of a set of beads connected together by massless segments of fixed length. Next, 
we compute the generating functional in the limit in which the chain becomes continuous. 
In order to do that, it is necessary to perform a path integral that is complicated by the 
presence of a functional Dirac delta function related to the constraints. The path integration 
is simplified by the introduction of suitable auxiliary fields. At the end, a nonperturbative 
expression of the generating functional is obtained. Finally, a similar method is exploited to 
compute the generating functional of the correlation functions of a two dimensional nonlinear 
sigma model with 0(3) group of symmetry. 



CONSTRAINED STOCHASTIC PROCESSES 

Let {R} = Ri{t), . . . , RN{t) denote the radius vectors of a set of particles in three 
dimensions. The positions of the particles are constrained by the conditions: 

C„({H}) = a= 1,...,M<3A^ (1) 

Their motion is subjected to the forces i^j = / + i/j for z = 1, . . . , A. Apart from a rescaling 
by a dimensional constant, namely the inverse of the friction coefficient, the = f^{{R}) 
are external forces, while the i^j's represent random forces with gaussian distribution, whose 
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properties are determined by the following correlation functions: 

{u,it)) = (2) 
{i^iAiHbit')) = 2D5,jSab5it-t') (3) 

Here a,b = 1,2,3 label the spatial components and D characterizes the width of the noise 
distribution. Assuming that the particles are fluctuating in a viscous solution, their motion 
will be described by the overdamped stochastic equation: 

R^ = fi + ly^ Z = l,...,iV (4) 

to which one has to add the conditions ([H). Following a standard notation, in Eq. (jll) it has 
been put Ri{t) = Before passing from Eqs. ([1]) and (jl]) to a path integral formulation 

of the generating functional of the correlation functions 

Cii,...,j„,ai,...,a„ (^1) • • • ) ^n) — (-Rii ,ai (^1 ) ; • • • j -Ri„,a,i (^n) ) (5) 



we briefly present a strategy explained in 411] in order to construct the generating functional 



in the simple case of the Langevin equation of a particle moving along the x— axis: 

mx + TJX — f = ly (6) 
with rj being the friction coefficient. Due to the inertial term, Eq. ([6]) is of second order in 



time. The idea of [4l[ is to transform it in a system of two equations of the first order by 
introducing the new variable 

V = X (7) 

In the enlarged two dimensional space {x,v), Eq. ((Tj) looks like a constraint which can 
be imposed a la Parisi-Wu by adding a fictitious random force u with a gaussian noise 
distribution characterized by the standard deviation D. In this way, one obtains from 
Eq. ([6]) two overdamped Langevin equations: 

mil + rjv — f = V (8) 
X — V = y (9) 

Starting from the above system of equations, the construction of the generating functional 
of the correlation functions is straightforward. Unfortunately, this method for fixing the 
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constraints cannot be applied to our case without changes because of two relevant differ- 
ences. First, the constraint ([7]) is imposed on a new variable, the velocity v, that was not 
present in the original Langevin equation ([6]). For this reason, there is since the beginning a 
clear distinction between the physical and redundant degrees of freedom. In principle, it is 
possible to extract also from Eqs. ([T]) and (jlj) a set of physical variables ui, . . . , u^n-m such 
that Ca(-Ri(Mi, . . . , U3n-m), ■ ■ ■ , Rn{ui, . . . , MsTv-Af)) = for « = 1, . . . , M. The velocities 
f 1 , . . . , vm could then be identified with the spurious degrees of freedom living in the sub- 
space of M.^^ which is orthogonal to that spanned by the m's. The major difficulty of this 
program is that, in order to determine the variables ui, . . . ,usn-m, it is necessary to solve 
the relations ([T]). This task is in general not feasible. On the second hand, Eq. ([71) is already 
in the form of an overdamped Langevin equation apart from the presence of the noise, while 
Eqs. ([1]) are not. To eliminate these differences, we introduce the new degrees of freedom 
and new random forces 9^, a = 1, . . . , M, characterized by the correlation functions: 

(/>a(t)) = (10) 

{i>a{t)Mt')) = ^DJ^pSit-t') (11) 
At this point, we replace Eqs. ([I]) and (jlj) with the following set of Langevin equations: 



Vaia + Ca{{R}) = 



(12) 

(13) 



In the equilibrium limit — )■ and by requiring that the friction constants rj^ vanish 
identically, the Langevin equations ( fT3ll reduce to the original constraints (JT]). The main 
advantage of having enlarged in this way the space of variables is that now it is possible to 
apply the standard procedure for constructing the expression of the generating functional. 
First of all, we introduce two sets of external currents Ji, i = 1, . . . , N and ja, a = I, ■ ■ ■ , M. 
The generating functional may be written as follows: 



N M 



TTTT lim lim / VuiVv^ 

1=1 a=l 







exp < 


{-IX 







N 



M 



i=l 



a=l 




X exp 



dt 



M 



N 



(14) 



_a=l i=2 

In the above equation the quantity tf denotes the upper limit of the time interval in which 
the evolution of the system is followed. We stress the fact that the current term for the 
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radius vectors i dt Yl!i=2 ' i'^i ~ Ji-i)-, with i denoting the imaginary unit, contains the 
discrete derivative of the external currents Ji and not the current itself. This definition of 
the current term is a just matter of future convenience and does not reduce the generality of 
our treatment. Let us also note that in Eq. the ^^'s and the iij's represent respectively 
the solutions of Eqs. (|T3l) and ( JT2l) and thus they depend on the noises Va and Ui. Moreover, 
the Riitys satisfy the time boundary conditions: 



Ri{0) — i?o,i 



R^{tJ) = R 



(15) 



Here the Ro/s denote the fixed initial conformations of the chain, while the Rf/s describe 
the conformation at the final time t = tj. Next, we perform in Eq. f|T^ the change of 



variables Ui , i>Q, 



Z[Jj] 



X exp 



WW lim lim / VRiV^a 



N 



dt 



M 

E 

a=l 



{:nic.+Cc({R})f 







j>exp < 









Jac exp \ — 

i=l 

V 

X] ^^^^ + i'Y^Ri- {Ji - J 



AD 



M 



a=l 



(16) 



where Jac is the determinant of the transformation. In principle, Jac is the determinant of 
a complicated block matrix: 



Jac = det 



(17) 



The determinant in the right hand side of Eq. f[T7j) can however be simplified by noticing 
that Eqs. (fT2|) establish relationships only between the i^j's and the R/s. As a consequence. 



it turns out that 



Sip 



for i 



,N, a 



1,2,3 and (3 



Jac reduces to a product of two determinants: 

Jac = det 



det 



M. For this reason. 



(18) 



A simple way to prove the above identity is to eliminate first the i/j's in the functional integral 
( |T^ using Eqs. ( |T2l) . Clearly, from those equations it turns out that Vi{R) = Ri — /j. The 



jacobian of this transformation is det |^ . After substituting everywhere in ( HM the t/j's 
with their new expression in terms of the -Rj's, we may eliminate also the auxiliary noises 
exploiting Eqs. f lT3|) . As a result, the Ua become functions of the new variables Ri^s and 



^Q,'s. The jacobian determinant of this second transformation is det 



Sip 



. This completes 
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our proof. Alternatively, Eq. (fT8|) may be easily checked by noticing that in the calculation 



of the full determinant in the right hand side of Eq. ( 1T7|) each element of the matrix 
necessarily multiplied by an element of the matrix which is equal to zero. 
Explicitly, the determinants in Eq. ( |T8l) read as follows: 



SR. 



IS 



det 



det 



with fi^ait) = /i,a(^i(i^), . . . , i?Ar(t)) and 



SijSab 



dt dRj^bit) 



5it - t'] 



det 



Va-^S{t-t') 



(19) 



(20) 



The only relevant determinant is that in Eq. (|T9|) . Determinants of this type can be evaluated 
in closed form. We give here only the result of the calculation referrin g th e interested reader 
for all details to standard books in this subject like for instance Ref. 10| : 



det 



''"dt dR,^,{t) 



6{t - t') 







= exp 


Hi 



N 3 



(21) 



Putting the expression of the jacobian just computed in Eq. (fT6l) and taking the limit rja — )■ 0, 
we obtain: 



Z[J] 



VR, 

1 JRi{0)=Ro,i 

t, N 3 



N 



exp 



—y 

AD ^ 



i=l 



X exp 

M 







^f^At) 



^^EE 







exp < 









N 



J i-l] 



1=2 



X 



JJ^ lim 



exp 



M 



a=l " 



dtClUR}) 



(22) 



The ^cf's decouple from the other degrees of freedom in the limit ?7q — )■ and have been 
already integrated out in Eq. ( l22l) . It remains to perform the limit Da — ?■ 0. To this purpose 
it is possible to apply the formula [3ll, |33[ : 

1 r*/ f2/ 



lim e-^ioV^{{K}) = 5(j({m)) 



(23) 



As a result, Eq. (1221) becomes: 
Z[J] - 



pR,{f)=Rf,, 



X exp 



7 Jiii(0)=Ho,i 



exp 



4D 



4 



g/..a(t) 

^.ttdRUt) 







exp < 









N 



i=2 



X 



M 



a=l 



(24) 
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This is the desired expression of the generating functional Z\J\. In writing the above formula 
we have have assumed that the forces j ^ are conservative, so that it is possible to use the 
identity fl^^ dt{Ri - f^f = j'l'' {R^ + f^) which is valid up to a term which is a total 
derivative. 



THE EXAMPLE OF AN INEXTENSIBLE CHAIN 



As an example of the method for imposing the constraints explained in the previous 
Section, we consider here the generating functional Z* '^ [J] describing the dynamics of an 
inextensible chain fluctuating in a viscous media. The chain is a discrete mechanical system 
consisting of beads of mass m connected together by — 1 massless segments of fixed 
length a. The total length of the chain is L = Na and its total mass is M = Nm. For this 
particular example the constraints ([1]) take the form: 

I-Rq, — -Rfj_ip — CL^ 



C^{{R}) 







a 



2,...,N 



(25) 



The direct application of Eq. ( I24l) to the inextensible chain gives as a result the generating 
functional: 



Z'-^-[J] 



N 



Riif) — Rf,i 



n / 



1 



N 



exp 



i=l 



X exp 

N 



I ftf ^3 







exp < 









N 



i=2 



X 



R. 



■a—l 



a=2 



(26) 



In the following all external forces will be switched off, i. e. /j = 0. We will now consider 
the continuous limit of Z* '^ [J]: 



a ^ 



^ 00 



Na 



const. 



(27) 



This limit can be performed using the prescription of Ref. 
function 

Z\J] = lim Z'-^-lJ 

a^O,N-^oo,Na=L 

of what has been called the generalized nonlinear sigma model in Ref. 3l[ |: 



3l|. The result is the partition 

(28) 



Z\J] 



R{tf,s)=Rf{s) 



R{0,s)=Ro(s) 



Pi2(t,s)e"^°''*^o'K^*''^+''^*''^-''M<5 {\R'it,s)\^ - 1) (29) 
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In the above equation Ro{s) and Rf{s) represent the continuous version of the discrete 
initial and final conformations of the chain given in Eq. f[T5|) . Moreover, s is the arc- length 
on the chain and R'{t,s) = ^^^i^^i. Finally, 



ds 

1 M 



(30) 



2kBTT 2L 

with kB being the Boltzmann constant. The parameter c contains physical constants like 
the temperature of the viscous medium T and the relaxation time of the infinitesimal beads 
r. In the frame of the duality between statistical physics and quantum mechanics, the 
quantity 2kBTT is the analog of the Planck constant [42]. It describes the uncertainties in 
determining the positions and the momenta of the beads because of the collisions with the 
molecules of the surrounding medium. 

Up to now, the generating functional Z[J] has been computed in the semi classical ap- 



proximation 311] or it has been linearized using a variational method 33|]. In the rest of 
this Section an attempt to go beyond these approximations will be presented. In order to 
avoid as much as possible complications with the boundary conditions, we will suppose that 
the chain is ring-shaped. Accordingly, throughout the rest of this Section all the fields will 
satisfy periodic boundary conditions of the type: 

R{t,s + L) = R{t,s) (31) 



The functional Dirac delta function imposing the constraints in Eq. fl29|) may be treated by 
introducing the new field r{t, s) and rewriting Z[J] as follows: 

Z[J]= / VR{t,s) ^^^^^^^^^f, dtf,[cR(t,s)+R(t,s).J(t,s)) 

J R{0,s)=Ro{s) J r{0,s)=R'f^{s) 

X 6{r{t,s)-R'{t,s))6{r\t,s)-l) (32) 

When t = and t = tj, r(t, s) obeys respectively the boundary conditions: 

r{0,s) = R[{s) and r{tf,s) = R'f{s) (33) 

At this point we express the first delta function appearing in the right hand side of Eq. (!32|) 
using the Fourier representation: 

S(r -R) = j Vke'^o^ ^'^^ dsk-{r-R') ^34^ 
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The new Fourier variable k = k{t, s) must be free to vary in the interval (—00, +00) over 
the whole domain [0, L] x (0, t/) of s and t in order to impose the relation r{t, s) = R'{t, s). 
When t = and t = tf, instead, the values of r(0, s) and r{tf,s) have been already fixed in 
the generating functional of Eq. ( l32l) to be those of Eq. ( l33l) . For this reason, at the instants 
t = and t = tf we may chose for k{t, s) Dirichlet boundary conditions, i. e. 

k{0,s) = k{tf,s) = (35) 

Applying Eq. ( IMIl in Eq. ( l32ll . it is possible to split 2[J] into a path integral over k{t, s) of 
the product of two different contributions: 

Z[J] = j VkZi[k,J]Z2[k] (36) 

where 

rR{tf,s)=Rf{s) ( r^f r ■ 2 l^ 

Zi[k,J]= / VRexp\- / dt ds cR + i{J' + k') ■ R \ (37) 

J R{0,s)=Ro{s) { Jo Jo ^ ^ J 

and 

Z2[k] = / Pre+'^o "^^^O dsk-r^^^2 _ (3g^ 

Jr{0,s)=R'o{s) 

As we see, the integration over R{t, s) is completely disentangled from that over r{t, s), but 
both fields interact with k{t, s). 

First, the functional Zi[k,J] will be computed. To this purpose, it is convenient to 
perform the change of variables: 

R= Rci + SR (39) 

Here Rd is a solution of the free classical equation of motion Rd = and satisfies the same 
boundary conditions as R. Consequently, the fluctuation 6R obeys the Dirichlet boundary 
conditions: 

5R{tf,s) = 6R{0,s) = (40) 

Moreover, both Rd{t, s) and 6R are periodic in s. After a gaussian integration over 6R, we 
obtain: 

J] = e-^-'C^'-^^exp/-^ dtdt' [ dsG{t,t'){J'{t,s) + k\t,s))-{J'it',s) + k'{t',s)) 
I 4 Jo Jo 

(41) 



where 

Sd{k, J) = dt ds 
Jo Jo 



cR,, + liJ' + k') ■ R^ 



(42) 
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while G{t, t') is the Green function satisfying the equation: 



and Dirichlet boundary conditions. Exphcitly: 



^c'^-^it-,, (43, 



2ctf 



(t + f) \t-t'\ 



(t + f) \t-t'\ 



(44) 



2 2 

Next, we consider the functional Z2[k]. This quantity closely resembles the Fourier trans- 
form of the probability distribution of a freely jointed chain appearing in the statistical 
mechanics of polymers. It may be computed exactly in a similar way using a discretiza- 
tion procedure. We have already seen how to pass from the continuous interval [0, L] to its 
discrete approximation in the case of the arc-length s, see e. g. Eq. (!27l) . Analogously, we 
discretize the time interval [0,t/] by replacing it with a lattice of Ai sites and link length b. 
The continuous limit is performed by requiring that: 

b^O M^+oo Mb = tf (45) 

Ignoring for the moment the nontrivial time boundary conditions of r(t, s), a straightforward 
calculation shows that Z2[k] is given by: 

N M 



Z.\k]= lim TTTT!!!^HfcJ (^g) 

N^^c^ Nn=L J- J- J- J- ab\k„rn\ 



a->0,Af^+oo,Afa=L 



Here knm is a shorthand notation for k{tm,Sn)- Apart from the boundary conditions, the 
problem with Eq. fH6|) is the evaluation of the continuous limit of its right hand side. This 



is a nontrivial task. The naive prescription used in the statistical mechanics of polymers 

sin(ab|fc I) ^ ^ _ ^^^|^^^|)2 ^ g-(a6|fc_|)^ (47) 

fcnm I 

fails to provide the correct result. For that reason, in order to evaluate Z2[k], we expand 
the exponential e'-^o'^'^*/o '^■^'"'^ in power series: 

^2|^] = ES^^2. (48) 

(T=0 



where 

ha= Vr{t,s)6{r\t,s)-1)T\ dtj dsika,{ti, si)ra,{ti, si) (49) 

Jr{0,s)=R'g{s) i^-^ Jo Jo 
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In the above equation the sum over repeated spatial indexes Oi, . . . , 020- = 1, 2, 3 is under- 
stood. Due to the presence of the functional delta function, it is difficult to apply the Wick 
theorem in order to integrate over the field r{t,s). However, a strategy to evaluate l2cr 
becomes clear after the discretization of both variables t and s. Calling l2a{N, A4, a, b} the 
discretized version of l2a, we obtain: 

N M-l ^ 2a / N M-l 



%imi T"ai,nimi 



n=l m=2 1=1 \n;=l m;=2 

TV 

JJ / dr^idrnMA^^M (50) 



N 

X 

n=l " 

with the boundary conditions at the initial and final times fixed by Dirac delta functions 
inside the factors Ai and Am- 

N / , ^\ N 

^ iM 



^Ro,n- Ro,n-l)\ A -Tf X I ^ - Rf,n-l, 

a 

(51) 



n=2 ^ a / 



Let us note that, due to the choice of Dirichlet boundary condition for k{t, s), the boundary 
values of this field are all zero, i. e. = knM = for n = 1, . . . ,N. For this reason, the 
sums over the indices mi in Eq. fl50|) have been restricted in the range 2 < mi < Ai — 1 for 
1 = 1,..., 2a. As a consequence, the variables at the boundary r„i and VnM in Eq. ( I50ll are 
not present in the product nf=i (Ei^i=i « E^^=2 ^ Kunirmrai^mm^- They enter only in the 
terms Ai and Am and are not mixed with the other degrees of freedom to be integrated. 
Thanks to this fact, the integrals over r„i and r„_A4 for n = 1, . . . , N can be easily factored 
out from the expression of l2a{N, JH, a,b}. The calculation of l2a{N, JH, a,b} requires the 
integration over the remaining 3A^(A^ — 2) variables r„m for 1 < n < and 2 < m < — 1 
with the measure drnm^ifnm ~ The integrand consists in a sum of terms in which 2a 
components ra^^mmi of the vectors ^'njm, with various combinations of the indexes ai,mi,ni 
are multiplied together. Each of such components may appear inside the product with 
powers of any order comprised between and 2a. Clearly, only the integration over terms in 
which the indexes ai of the components and the discrete indexes ni,mi are present in even 
combinations will not vanish identically. Basing ourselves on these observation, it is possible 
to establish the following rules in order to evaluate /2o-{A, A^, a, b}: 

1. Rearrange the terms appearing in l2a{N,J^, a, b} in such a way that the products of 
the 2a components rai^mmi^ will be grouped into products of 2/— plets for / = 1, . . . , cr 
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characterized by the fact that inside each multiplet all the components have the same 
indices ni and m^. 

2. Eliminate the spatial indices ai by using the formula 

3 

^ ^ ^ai,nimi'^ai,nimi | ^n;m; 1 1 '''njm; | COS ^j^^mj (^^) 



3. Use the fact that, due to the delta functions S{rl^^ — 1) it is possible to put |r„jm; | = 1 
and integrate over the angles Ommi using the formula: 
r 1 

/ d6mn sin 6nm COS^* 9nm = 2 i = 0, 1, . . . (53) 

Jo 2t + l 

Let Ki be the number of 2/— plets mentioned in [H Clearly, the possible values of Ki are 
limited by the condition Yl'i=i '^^^i = 2(T. This implies for instance that Ki = 1 for I > a. 
Moreover, the number of combinations for dividing 2a objects in a number Ki of pairs, K2 
of 4— plets etc. is given by: 

2(t! 

^ " 2Ki!4K2!---2(a-l)K._i! ^^'^^ 
Following steps [H and[2l of the above prescription, we may write: 

, a N M 

hAN,M,a,b} = y ■ , -TT y a^' V 6^' 

Ki,...,Ka «=1 x,=l m.i,---,mi,K,=^ 

Er=i 2iK,=2a 

^1^ lA/'U |2_1^|^ |2/ 1^ 

l.KilTT-l.Ki I 



^ l^ni,i»Tii,i I • • • l^n;^K;"i;^K; I / '^|^nr?i|'^(|''^n,m| ~ 1) I '''ni^iJTii^i | • • • l^n;^^'; 

n,m=l 

/'27r /"TT 

X / d(pnm / C^^'nm slu 6'„m COS^' K,imi^l ■ ■ ■ COS^' 6ni^j^^m,,K,'^l'^M (55) 

It is now possible to perform the remaining integrations according to prescription [3l The 
result is: 

2a\ 



Kl,...,Kcr 



N M 

X 
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FIG. 1: This figure shows the Feynman diagrams corresponding to the terms J^^ dt ds\k{t, s)p, 



dt ds\k{t,s)\'^, ff^^ dtf^ ds\k{t,s)\^,. . . appearing in Eq. ([58]) . 
where 

/ de 

n,m=l -^0 -^0 

is the discrete version of Z2[k = 0]. The continuous hmit of l2cr{N, Ai, a,b} starting from 
Eq. gives: 

r , 2o"! 



^ ^ 2Ki!4K2!---2(a-l)K._i! 

Kl,...,Kcr 



X 



n (2/ + ^^''1 dsi^i\Kki,si^i)f ■■■ ^^''^'/o dsi^KAKkKv si,Ki)?^ (58) 

with 

ZsfOl = hm ZafO; A^,7W,a,6] (59) 

b^O,M^+oo,Mb=tf 

Z2[0] contains the dependence on the boundary values of the fields R{t,s). Finally, sub- 
stituting Eq. ( 158]) in Eq. ( 148|) . we obtain the explicit expression of Z2[k]. Having derived 
the functionals Zi[k, J] and Z2[k], we may compute the full generating functional Z[J] of 
Eq. dMl): 



-oo 



^ 2a\ ^ 2Ki!4K2!---2(a- l)AV-i! (2/ + 1)^^' 



X / / dsi^i\k{ti^i,si^i)f ■ ■ ■ I dti^Ki [ ds^Ki\H'ti,Ki, si,Ki)f (60) 

^0 Jo Jo Jo 

and perform the remaining integrations over the fields k{t, s). Let us note that the action 
in Eq. (!60|) is gaussian. The classical sector present in the term Sci{k, J) of Eq. (H2|) intro- 
duces a coupling of k{t, s) with an external current which is proportional to the classical 
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conformation Rci{t, s). It is actually possible to rewrite the generating functional 2[J] in a 
much more compact form by computing the following functional: 

f^f ■ 2 

Z[J,fi] = exp{— / dt dscR^i{t, s)} 
Jo Jo 

X exp |-* ^ ^ dtdt' dsG{t, t'){J'{t, s) + k'{t, s)) ■ {J'{t', s) + k'{t', s))| 
X exp I ~i I dt I dsfj.{t, s) ■ k(t, s) I (61) 



where fi{t, s) is the external current that is needed to generate the correlation functions of 
the field k{t, s). Z[J, fi] is the double generating functional of the correlation functions of 
the physical fields R{t, s) and of the auxiliary fields k{t, s). Knowing Z[J, fi], the generating 
functional Z[J] may be expanded as follows: 

+ 00 \ 1 

^[•^] = 1] n (2lKi)\ (21 + 

cr=0 Ki,...,K^ 1=1 ^ ^ ^ 

j:i^^2lKi=2a 

\6fi{ti^i,sm) ■ 6fi{ti^i,si^i) ) \6iJ.{ti^K,,si^K,) ■ S^j,{t^Ki,si^Ki) ) ^'"^'^l 1/^=0*^^^^ 

Z[J, fi] can be computed in closed form. Its expression is given below: 

Z[J,fi] = Z2[0]exp|— y dt J ds (^cit^^i + iJ' ■ Rci — tfJ, ■ J 

X exp I - [ ^ dtdt' [ dsds'g{t, s; t' , s')il{t, s) ■ jl{t', s') } (63) 

where 

a(t R- f' r'\ = 

2 ^ (27m) 2 L L 

n=l ^ ' 



, . cd'^bit-t')^ L2 . 2'Kns . 2W 



and 

fi = fi + Rd (65) 

It is straightforward to realize that at each order in the index a, the series in Eq. (162|) contains 
products of vertices with 21 external legs like those shown in Fig. [1] for / = 0, 1, 2, . . . , cr. 
The legs of the 2/— vertices are contracted in all possible ways in Eqs. fl60p or fl62p using the 
propagator flM|) . Convergence is granted for small values of the constant c. It is important 
to stress that the term 22[0], which contains a nontrivial dependence on the boundary 
conformations of R{t, s), appears in Eq. (!63|) as an overall factor that can be easily eliminated 
by choosing the normalization of the correlation functions. 
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A BRIEF DIGRESSION ON THE NONLINEAR SIGMA MODEL 

In this Section we will apply the method illustrated in the case of the inextensible chain 
to a two dimensional nonlinear sigma model with action: 

Ssm[J] = [ d^x [| {d,ct>f + Jc/)] (66) 

where di = x = {xi, X2) and M is a two dimensional manifold with Euclidean signature. 
To avoid complications with the choice of boundary conditions, M is chosen to be a torus, 
so that periodic boundary conditions should be implemented for xi and X2- The vector field 
4){x) is subjected to the constraint: 

0^(x) = 1 (67) 

The generating functional of the nonlinear sigma model in path integral form looks very 
similar to that of the inextensible chain: 

Zsm[J] = J Vct>e-'-^'^{cf>^ - 1) (68) 

Introducing the new fields r{x) and fc(x), the integration over the physical field (p can be 
performed without worrying about the functional Dirac delta function 5{(f)^ — 1)- The result 
of this integration is: 

Zsm[J] = J Vk{x)e~-^^M'i'-^'yGi^'y)(J(-)+H-))<Ji-)+H^))z^[k] (69) 

Apart from the different topology of the underlying manifold M, Z2[k] is the same as the 
functional Z2[k] given in Eq. 



Z2[k] = j Vr{x)e'^MH^yri^)s^r\x) - 1) (70) 

One may thus evaluate this functional exactly as in the previous Section. The result is an 
expansion of Z2[k] at all orders in powers of k{x) which converges for small values of the 
coupling constant g. The most relevant quantity to be computed in order to derive term by 
term at all orders the contributions to this expansion is the extended generating functional 
Zsm[J, analogous to that of Eq. (!6T]) . After a long but straightforward calculation, one 
obtains apart from a trivial overall factor the following expression of Zsm[J, fj]'- 
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where 



g{x,y) = '-^j\x,y) (72) 



and = #7 + 



CONCLUDING REMARKS 

In this work a recipe has been presented for constructing the generating functional of 
constrained stochastic systems with an arbitrary number of degrees of freedom and general 
constraints. The main idea is that the variables to be constrained are allowed to fluc- 
tuate thanks to the addition of the auxiliary noises z/q, and of the fictitious variables ^q, 
a = 1, . . . , M. In this way, the conditions Ca{{R}) in Eq. ( !T3|) are not fixed exactly, but are 
subjected to stochastic fluctuations. The original constraints ([1]) are recovered in the "rigid" 
limit r]a, Da 0. The set of equations ( 1T2|) and ( 1T3|) which describes the dynamics of both 
physical and fictitious degrees of freedom is formally a system of overdamped Langevin equa- 
tions. Starting from a system of that kind it is possible to derive a path integral expression 
of the related generating functional by using the standard techniques valid for stochastic 
systems without constraints. In the rigid limit r]^, = 0, the generating functional Z[J] 
is given in Eq. flM|) . Remarkably, with this procedure the jacobian determinant appearing 
after the change of variables Ui,i)a — > Ri,^a remains relatively simple. For comparison. 



the standard treatment of constrained path integrals, see for instance [36|, produces com- 
plicated jacobian determinants of block matrices. The path integral expressions obtained in 
this way are very useful at a theoretical level, because they allow for example to check if 
the Parisi-Wu quantization scheme delivers in the equilibrium limit the desired generating 
functional of the theory to be quantized js, 9|, but are very cumbersome in concrete calcula- 
tions of physical observables. The method discussed in this work produces instead relatively 
simple generating functionals. In the case of the inextensible chain in the absence of external 
interactions, we obtain for instance the GNLcxM of Ref. 
of several measurable quantities. For example, the dynamical form factor of the chain has 
been evaluated in 4J]. Always in 4^ it has been estimated how the fluctuations of the 



3l| . which allows the computation 



distance between two arbitrary points on the chain are influenced by physical parameters 
like the length of the chain and the relaxation time r. The novelty of the present approach 
with respect to [31] is the possibility to add also the external forces f^. 
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Of course, the expressions of the generating functionals obtained here are still complicated 
due to the presence of the Dirac delta functions that are needed to impose the constraints. 
For this reason, in this work a technique is developed for computing the the generating 
functional of two dimensional field theories in the presence of a Dirac delta functions like 
that appearing in the GNLaM of Eq. (l29l) . Thanks to the introduction of the auxiliary field 
r{t, s), the difficulties related to the delta function are confined to the computation of the 
path integral Z2[k] of Eq. f p8|) . Apart from the boundaries at t = and t = tf, Eq. f l46|) 
shows that Z2[k] formally consists of the infinite product over the time index m of terms 
that are nothing but the Fourier transform of the probability function of a freely jointed 
chain in statistical mechanics. The functional Z2[k] has been evaluated here in the form of 
a series converging for small values of the parameter c given of Eq. ( 130|) . Despite the fact 
that Eq. fl62|) looks like a perturbative expansion, all terms entering in the expansion can be 
computed in closed form. Indeed, the main ingredient in Eq. ( 162|) is the extended generating 
functional Z[J,fi\, which generates both correlation functions of the physical fields R{t,s) 
and of the auxiliary fields fc(t, s). The exact expression of Z[J, /j] is given by Eq.( |63l) . 

Finally, the techniques presented in this work can be extended to other systems with con- 
straints. For instance, we have seen there are many similarities between the two dimensional 
nonlinear sigma model and the GNLcrM. Despite the initial difference of the constraints, af- 
ter the introduction of auxiliary fields the generating functional of both theories may be 
computed in the same way. In particular, the extended generating functional ZsmiJ, fj] of 
the correlation functions of the physical fields </>(x) and the auxiliary fields k{x) for the 
nonlinear sigma model has been exactly derived in Eq. ( |7T]) . We hope that in the future 
it will be possible to apply these methods also to the statistical mechanics of polymer in- 



teractions in the presence of topological entanglement |45l-l47l| . In statistical mechanics, in 
fact, the probability function of a system of long and flexible polymer chains obeys pseudo- 
Schrodinger equations that, for certain aspects, are similar to those of Fokker-Planck. A 
nice physical application could be the study of the influence of topological constraints on 



the interactions in colloid-polymer mixtures j48|. 



49|. 
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